Semiclassical Prediction for Shot Noise in Chaotic Cavities 
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We show that in clean chaotic cavities the power of shot noise takes a universal form. Our 
predictions go beyond previous results from random-matrix theory, in covering the experimentally 
relevant case of few channels. Following a semiclassical approach we evaluate the contributions of 
quadruplets of classical trajectories to shot noise. Our approach can be extended to a variety of 
transport phenomena as illustrated for the crossover between symmetry classes in the presence of a 
weak magnetic field. 

PACS numbers: 73.23.-b, 72.20.My, 72.15.Rn, 05.45.Mt, 03.65.Sq 



X 



Ballistic chaotic cavities have universal transport prop- 
erties, just as do disordered conductors. The explanation 
of such universahty cannot rely on any disorder average 
but must make do with chaos in an individual clean cav- 
ity. We shall present here the semiclassical explanation 
of shot noise, relating the quantum properties of chaotic 
cavities to the interference between contributions of mu- 
tually close classical trajectories. Similar methods have 
recently been used for explaining universal spectral fluc- 
tuations of chaotic quantum systems 0,0, and to calcu- 
late the universal mean conductance in |3|, 3 . 

Following Landauer and Biittiker IHHi '^^ treat trans- 
port as scattering between two leads attached to the cav- 
ity. One lead is assumed to support A'^i ingoing channels 
and the second one N2 outgoing channels. In contrast 
to the random-matrix treatment of Q and work on 
quantum graphs in ^ , our results cover all orders in the 
inverse number of channels, = A^i -I- N2, and thus ap- 

fto the experimentally relevant case of few channels 
Previously unknown and surprisingly simple expres- 
sions for the shot noise arise, both with and without time 
reversal invariance (see Eq. Ullfl below) . 

The transition amplitudes between ingoing channels 
ai and outgoing channels 02 define an A'^i x N2 matrix 
t — {taia^}- That matrix determines the power of shot 



noise as 



P = (tr(ttt _ ttht^)), in units 



2e^|V| 



depending 



on the voltage V; for us, (. . .) denotes an average over 
a small energy interval. Previous work had involved av- 
erages over ensembles of matrices t and obtained 



P = — — - -\- 

N3 



1 



N,N2iN,-N2)^ 
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here /3 — 1 refers to the so-called "orthogonal case" of 
time-reversal invariant dynamics; if a magnetic field is 
applied to break time- reversal invariance ( "unitary case" , 
/3 = 2), the second ("weak localization") term disap- 
pears. Higher orders in -i- are as yet unknown. 

In the semiclassical limit, each transition amplitude 
taia2 is given by a sum over trajectories a leading from 
an ingoing channel ai to an outgoing channel 02, iaia2 ~ 
Ea(ai-.a2) ^^e''^"/'* It Can be shown that the ab- 
solute value of the initial angle of the relevant trajectories 



(i.e. the angle enclosed between the initial piece and the 
direction of the lead) is dictated by the ingoing channel, 
whereas the final angle is determined by the outgoing 
channel. The contribution of each trajectory depends on 
the Heisenberg time Th = (2^^)/'-! ' "^ifi^ ^ f^e volume 
of the energy shell and / the number of freedoms. The 
factor Aa is determined by the stability of the trajectory, 
and the phase is proportional to the classical action ■ 
With the transition amplitudes thus semiclassically ap- 
proximated, the quadratic term (tr(tt^)) becomes a dou- 
ble sum over trajectories. That double sum, which ac- 
tually is the mean conductance, was evaluated in 4] as 
N -1+2/13 • "^^'^ quartic contribution to shot noise turns 
into a sum over quadruplets of trajectories 



(tr(tt^tt^)) = 



E 
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here ai, ci — 1, . . . , iVi and 02, C2 — 1, . . . , N2 represent 
ingoing and outgoing channels, connected by the trajec- 
tories a, /3, 7, 6 like a (ai — > 02), /3 (ci — > 02), 7 (ci C2), 
6 (ai — > C2). The sum is dominated by quadruplets where 
the trajectories /3 and S have approximately the same 
cumulative action as a and 7, such that the action dif- 
ference AS" = Sa — Sp ~ Ss is of the order h. The 
contributions of other quadruplets interfere destructively. 

Diagonal contribution: The simplest quadruplets have 
either a — (3,j = S,0Ta — S, (3 — j 8J; their action 
difference vanishes. The first case has coinciding ingoing 
channels ai and Ci and contributes 



(tr(tt^tt^))o=,3 

7 = 5 
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The foregoing sum can be done using ergodicity. As 
shown in Q, summing all trajectories between two fixed 
channels amounts to integrating over the dwell time T, 



E 1^^ 
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where e ^« is the probabiHty for the trajectory to dwell stretch. The probability that both a and 7 stay inside 



in the cavity up to the time T, and ^ the rate of escape. 

To proceed with Eq. (|3Jl we invoke the Richter/Sieber 
sum rule twice and afterwards sum over all N1N2 
possible combinations of channels with ai — ci. Simi- 
larly, the case a = (5, /3 = 7 leads to N^N2 combina- 
tions with coinciding outgoing channels 02 = C2. Alto- 
gether, these so-called diagonal contributions sum up to 
_ . In the unitary case, they cancel 

with (tr(tt''')) such that shot noise must be entirely due 
to different quadruplets of trajectories. 

2-encounters: The first family of such quadruplets, de- 
picted in Fig. ^ was identified by Schanz, Puhlmann and 
Geisel for quantum graphs j^]- Here, the trajectories a 
and 7 approach each other in a "2-encounter" : a stretch 
of a comes so close in phase space to a stretch of 7 that 
the motion over the two stretches is mutually linearizable. 
The remaining parts of a and 7 will be called "loops". 
Assuming that all loops have non-vanishing length (En- 
counter stretches do not "stick out" into the leads), one 
finds two further trajectories /3 and (5, which practically 
coincide with a and 7 inside the loops but are differently 
connected in the encounter: The trajectory f] closely fol- 
lows the initial loop of a and the final loop of 7, whereas 
6 follows the initial loop of 7 and the final loop of a. Ob- 
viously, the cumulative action of (3 and S approximately 
coincides with the action of a and 7, with the action dif- 
ference exclusively determined by the encounter region. 




FIG. 1: Quadruplet of trajectories a, (3, 7, 5 differing by their 
connections inside a 2-encounter (in the box) . Initial and final 
points marked by channel indices ai , ci , 02 , C2 . The durations 
are tcnc for the encounter and ti, t2, tg, ti for the loops. 

Each encounter infiuences the survival probability. If 
a particle stays in the cavity along the first encounter 
stretch it cannot escape during the second stretch either, 
since the two stretches are close to each other. The tra- 
jectories a and 7 are thus exposed to the danger of get- 
ting lost only on the four loops (see Fig. ^) and on one 
encounter stretch. Denoting the duration of the latter 
stretch by tone wc can write the overall exposure time as 
T^exp = ^1+^2+^3 + ^4 + ^enc- That cxposure time is 
smaller than the cumulative duration Ta + oi a and 
7, by a second summand tcnc for the second encounter 



the cavity reads e ^« ""'^ , larger than the naive estimate 
Q~T-ii'^'^°'~^'^''\ In brief, encounters hinder escape Q. 

To describe the geometry of encounters, we consider a 
Poincare section V in the energy shell, through an arbi- 
trary point of a. If 7^ cuts through an encounter as in 
Fig. n it must intersect 7 in a point close to the refer- 
ence point on a. Assuming two freedoms we decompose 
the separation between both points into coinponents s, 
u along the stable and unstable manifolds Both s 
and u must be small, |u| < c, \s\ < c, with c some clas- 
sically small constant. The components s and u fix the 
action difference as AS* = su and the encounter duration 
as ienc = X I where A is the Lyapunov constant |2[ . 

Using ergodicity, we count the encounters within tra- 
jectory pairs. The probability density for 7 to pierce 
through P at a specified time with phase-space separa- 
tions s and u is uniform and given by the inverse of the 
volume of the energy shell Q. To capture all encounters, 
we integrate that density over (i) the time of piercing on 
7 and (ii) the time of the reference piercing on a. The 
probability for V to cut an encounter is proportional to 
the duration t^^, which we divide out to get the number 
of encounters [21 • Changing the integration variables to 
the loop durations ti and ^3 we get the weight function 



w{s, u) 



dti 

^0 



1 



; (5) 



riienc(s,1i) 

here the upper boundaries, depending on the dwell times 
Ta, Try of a and 7, make sure that the loop durations t2 
and ti remain positive. The density w{s, u) is normalized 
such that J dsduw{s,u)S(su — AS) is the number density 
of 2-encounters with fixed action difference AS". 

To account for all quadruplets of Fig.^ we do the sum 
over /3, (5 in (0) by integrating with the weight w(s,m), 



T^ 



a, 7 



approximating AjsAs « AaA^ j^. Now, similarly to 
@, we replace the sum over a, 7 by an integral over the 
dwell times or, equivalently, over the loop durations t2 
and t^. The integrand must be weighted with the prob- 

ability e ™p for both trajectories to remain inside 
the cavity. Regarding the sum over channels ai, 02, ci, C2 
one might expect a factor iV^A''|. However, when both 
the ingoing channels and the outgoing channels coincide, 
fli = ci, 02 = C2, each of the two dashed trajectories in 
Fig- n could be chosen as /3 or S, such that the resulting 
contributions are doubled. Including such combinations 
for a second time, we get the factor NiN2{NiN2 + 1). 
Altogether, we thus find 



(tr(tt^tt^)) 



2— cnc 



_ NiN2{NiN2 + l) 



J dtidt2dtzdti 



N 



[tl+i2+i3+t4 + to„c(s,«)]„i|2i 



^). (7) 
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FIG. 2: Families of trajectory quadruplets a,/?, 7, 5 responsible for the next-to-leading ("weak localization") contribution to 
shot noise (drawn similar as in Fig. but without the cavity). Each picture represents either four or two similar families of 
quadruplets. Arrows indicate the direction of motion inside the encounters (thick lines), and highlight loops which are traversed 
in opposite direction by (a, 7) and {/3,5). 



The integral factors into four independent integrals over 

the loop durations, Jp dtie ' = and one in- 

tegral over the separations s, u inside the encounter, 



as shown in 



/ ds duj 

0. Since all powers of Th mutually cancel, the following 
diagrammatic rule arises: Each loop gives rise to a fac- 
tor jj , and an encounter contributes a factor —N. The 
rule yields upon multiplying with the number of 

possible combinations of channels we get 



-{tr{ttht^))2-c 



N^N^jN^N^ + 1) 
N3 



(8) 



i.e., for Ni, N2 ^ 1 the leading shot-noise term in 

All orders: To go beyond the leading term, we have 
to account for all quadruplets of trajectories differing in 
arbitrarily many encounters, each involving arbitrarily 
many stretches; Fig. |5] shows a few examples. In the 
unitary case we must consider encounters where several 
stretches of either a or 7, or both, come close in phase 
space, whereas in presence of time- reversal invariance the 
stretches may also be nearly mutually time-reversed. 

The contributions of all families of quadruplets obey 
the above diagrammatic rule. To show this, we describe 
each Z-encounter (encounter of I stretches) by Z — 1 pairs 
of coordinates sj, uj, j = 1, . . . , Z — 1 Q measuring the 
separations of Z — 1 stretches from one reference stretch. 
These coordinates determine both the duration of each 
encounter stretch and its contribution ■^j'^j to the 

action difference. The analog of the density w{s, u) in (O 



obtains a factor 



from each Z-encounter. The 



resulting product must be integrated over the durations 
of all loops, with integration over the final loops of a 
and 7 coming into play through the summation over a 
and 7 like in (0 . Finally, the contribution of each family 
factors into "loop" and "encounter" integrals similar to 
those in |7J. After cancellation of all powers of Tff, the 
diagrammatic rule comes about, with a factor from 
each loop and a factor —N from each encounter. 

Again, we have to multiply the result with the number 
of possible combinations of channels. Two cases must be 
distinguished. First, let us consider trajectory quadru- 
plets as in Fig.|2t-c where similarly to Fig.^the partner 
trajectories /? and 5 connect the initial point of a to the 



final point of 7, and the initial point of 7 to the final point 
of a. Such quadruplets will be called x- quadruplets. For 
them, the channels oi, ci, 02, and C2 may be chosen ar- 
bitrarily and allow for NiN2{NiN2 + 1) combinations. 

In contrast. Fig. |3i depicts a quadruplet where the 
partner trajectories connect the initial point of a to the 
final point of a, and the initial point of 7 to the final point 
of 7 , similarly to the diagonal contribution. We speak of 
a d-quadruplet then. The partner trajectories now con- 
nect the leads as ai — > 02, ci — > C2. Since for shot noise 
we need partner trajectories /3 (ai ^ C2) , (5 (ci 02), d- 
quadruplets contribute only if either the two ingoing 
channels, or the two outgoing channels (and thus the 
corresponding angles of incidence) coincide. Like for the 
diagonal contribution, we thus obtain NIN2 + NiN^ = 
NN1N2 possible combinations. 

Our diagrammatic rules determine {tY{ttHt^)) as 



{ti{ttht^)) 



NiN2iNiN2 + 1) ^ X, 




(9) 



Towards explaining Xm and dm we denote the number of 
encounters in a quadruplet by V and the total number 
of the encounter stretches by L. Then a;,„ is the number 
of families of x-quadruplets with m — L — V and even 
V, minus the number of corresponding families with odd 
V; dm is the analogous number of d-families. We note 
that the contribution of each family is proportional to 
rather than j^, since there are two more loops 
than encounter stretches. 

The leading contribution to shot noise originates from 
the family of Fig. ^ it gives xi = —1. For the next- 
to-leading term, we have to consider x-quadruplets with 
two 2-encounters or one 3-encounter, contributing to 
X2 and depicted in Figs. 2a-c, and d-quadruplets which 
are related to a single 2-encounter and contribute to di 
(Fig. 131). All quadruplets in Fig. 2 involve mutually 
time-reversed loops and can exist in the orthogonal case 
only (thus 0:2 = di = in the unitary case). Note that 
if we interchange the two leads, or the pairs (0,7) and 
(/3, (5), or the trajectories a and 7, each family of quadru- 
plets will be either left topologically invariant or turned 



4 



into an equivalent family making the same contribution 
to the shot noise. Only one representative of each such 
"symmetry multiplet" is shown in Fig. El with the num- 
ber of equivalent families indicated by a multiplier. The 
sum of all contributions gives X2 = 4, c?i = —2, and 
(tr(tttttt)) ^N^^^+ 4^ - 2^ 
Together with {tr(«t)) = „ 
cover the second term in 

For higher orders in we must collect all families of 
trajectory quadruplets. We had previously established a 
method for counting families of pairs of periodic orbits 
differing in encounters, based on permutation theory ■ 
The families of orbit pairs thus obtained can be turned 
into the families of trajectory quadruplets needed now, 
simply by cutting each pair twice, inside loops. One can 
show that if one and one only of the loops cut is tra- 
versed in the opposite sense in the orbits of the pair, the 
resulting quadruplet is of type x, and otherwise of type 
d. Using this method we obtain 



(-1)" 



2 unitary 
(_l)™ 3'"-i orthogonal 



unitary 



(-1)"^^ orthogonal. 



(10) 



changing direction in {(3,5) relative to (a, 7). Our above 
diagrammatic rule is thus modified: Each loop chang- 
ing direction contributes with ^ oc ; each en- 
counter contributes —N{1 + ■Cm^ number of 
its stretches changing direction 

The leading contribution to shot noise, from quadru- 
plets as in Fig. ^ remains unaffected by the magnetic 
field. However, all quadruplets responsible for the next- 
to-leading term obtain a Lorentzian factor (In 
Figs.|21i,[2t, and|2Ji one loop is traversed in time-reversed 
sense and /i = 0, whereas in Fig. two loops are time- 
reversed and one encounter has ^ = 1.) We thus predict 



P = 



N3 



O 



N 



(12) 



in accordance with pi|l for the limits ^ and ^ ^ 00. 
If Ni — N2 = N/2 the second term in H12|) vanishes, and 
quadruplets involving more encounter stretches yield 



16 N 



1 1 + 8e + + + 
16(1 + 0^ 



(13) 



The proof, based on a recursion derived in will be 
given elsewere. Summing over m we get the shot noise 



N{N^-1) 

Ni{Ni + l)N2{N2 + l) 



unitary 
orthogonal , 



(11) 



valid to all orders in and thus also for a few channels. 

Weak magnetic field: In the presence of a weak mag- 
netic field B, the power of shot noise must interpolate 
between the orthogonal and unitary cases. A weak field 
increases the action of each trajectory 0, ^| by the line 
integral i J A • dq of the vector potential A. Since that 
increment changes sign under time reversal a net con- 
tribution survives from all loops and encounter stretches 



the extension to 0{^) will be given elsewhere. 

Outlook: The semiclassical approach to transport 
opens a large field of experimentally relevant applica- 
tions. For instance, we have checked the trajectory 
quadruplets employed here to also yield the universal 
conductance variance as well as the covariance of the 
conductance at two different energies, both within and in 
between the orthogonal and unitary symmetry classes. 

We are indebted to Dmitry Savin, Hans-Jiirgen 
Sommers (who have reproduced our prediction Hll|) 
in random-matrix theory), and Martin Zirnbauer for 
useful discussions and to the Sonderforschungsbereich 
SFB/TR12 of the Deutsche Forschungsgemeinschaft for 
financial support. 
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